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Abstract 

A new approach to distributed cooperative beamforming in relay networks with frequency selective 
fading is proposed. It is assumed that all the relay nodes are equipped with finite impulse response 
(FIR) filters and use a filter-and-forward (FF) strategy to compensate for the transmitter-to-relay and 
relay-to-destination channels. 

Three relevant half-duplex distributed beamforming problems are considered. The first problem 
amounts to minimizing the total relay transmitted power subject to the destination quality-of-service 
(QoS) constraint. In the second and third problems, the destination QoS is maximized subject to the 
total and individual relay transmitted power constraints, respectively. For the first and second problems, 
closed-form solutions are obtained, whereas the third problem is solved using convex optimization. 
The latter convex optimization technique can be also directly extended to the case when the individual 
and total power constraints should be jointly taken into account. Simulation results demonstrate that 
in the frequency selective fading case, the proposed FF approach provides substantial performance 
improvements as compared to the commonly used amplify-and-forward (AF) relay beamforming strategy. 
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I. Introduction 

Recently, cooperative wireless communication techniques gained much interest in the literature as they 
can exploit cooperative diversity without any need of having multiple antennas at each user [l]-[9]. In 
such user-cooperative schemes, different users share their communication resources to assist each other 
in transmitting the information throughout the network by means of relaying messages from the source 
to destination through multiple independent paths. 

Different relaying strategies have been proposed to achieve cooperative diversity. Two most popular 
relaying strategies are the amplify-and-forward (AF) and decode-and-forward (DF) approaches [1], [4], 
[6]-[ll]. In the AF scheme, relays simply retransmit properly scaled and phase-shifted versions of their 
received signals, while in the DF scheme the relay nodes decode and then re-encode their received 
messages prior to retransmitting them. Due to its low complexity, the AF relaying strategy is of especial 
interest [12]-[18]. 

When the channel state information (CSI) is not available at the relay nodes, distributed space-time 
coding can be used to obtain the cooperative diversity gain [10]-[13]. However, with available CSI, 
distributed network beamforming can provide better performance [14], [15]. 

Recently, several distributed AF beamforming techniques for relay networks with flat fading channels 
have been developed [14]-[18]. The approaches of [14]-[16] optimize the receiver quality-of-service 
(QoS) subject to the individual and/or total power constraints under the assumption that the instantaneous 
CSI is perfectly known at the destination or relay nodes. In these approaches, the QoS is measured in 
terms of the signal-to-noise ratio (SNR) at the receiver. In [17], the problem of AF relay beamforming 
is considered under the assumption that the second-order statistics of the source-to-relay and relay- 
to-destination channels are available. Based on the latter assumption, several distributed beamforming 
algorithms are developed in [17]. In the first technique of [17], the total relay transmit power is minimized 
subject to the receiver QoS constraint, whereas in the second approach of [17], the receiver QoS is 
maximized subject to the total or individual relay power constraints. In [18], the approach of [17] has 
been extended to multiple source-destination pairs. Recently, the problem of using an imperfect (e.g., 
quantized) CSI feedback in distributed beamforming has been also considered [19]. 

Although some extensions of distributed space time-coding techniques to the frequency selective fading 
case are known in the literature [20], the problem of distributed beamforming in frequency selective 
environments has not been addressed so far. In particular, all the techniques of [14]-[19] assume the 
transmitter-to-relay and relay-to-destination channels to be frequency flat. However, in practical scenarios 
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these channels are likely to be frequency selective. In the latter case, there is a significant amount of 
inter-symbol interference (ISI) which makes it difficult to directly extend the techniques of [14]-[19] to 
frequency selective fading channel scenarios. 

In this paper, we consider a relay network of one transmitter, one destination, and multiple relay 
nodes under the assumption of frequency selective, finite impulse response (FIR) transmitter-to-relay and 
relay-to-destination channels. To compensate for the effect of such channels, a new filter-and-forward 
(FF) relaying protocol is proposed. According to the FF strategy, all the relay nodes are equipped with 
finite impulse response (FIR) filters that are used to compensate for the transmitter-to-relay and relay- 
to-destination channels. 

Three relevant distributed beamforming problems are considered under the assumption that the instan- 
taneous CSI is available at the receiver or at the relay nodes. Similar to the techniques of [14]-[18], the 
receiver is assumed to use a perfect source-to-relay and relay-to-destination CSI to compute the relay 
weight coefficients and feed them back to the relay nodes using a low-rate receiver-to-relays feedback 
link. Alternatively, the relays can directly compute their weight coefficients, provided that the CSI is 
available at the relay nodes. 

Our first distributed beamforming problem amounts to minimizing the total relay transmitted power 
subject to the destination QoS constraint. As in the frequency selective case the major performance limiting 
factor is ISI rather than noise, the QoS is measured in terms of the receiver signal-to-interference-plus- 
noise ratio (SINR), in contrast to the techniques of [14]-[17] that use the SNR as a measure of QoS in 
the flat fading case. In our second and third problems, the destination QoS is maximized subject to the 
total and individual relay transmitted power constraints, respectively. For the first and second problems, 
closed-form solutions are obtained, whereas the third problem is solved using convex optimization. The 
latter convex optimization technique can be also directly extended to the case when the individual and 
total power constraints should be jointly taken into account. 

It is shown that in the flat fading case, the proposed FF network beamforming techniques reduce to 
the AF network beamformers of [14]-[17] which are particular cases of our techniques. 

In frequency selective channel scenarios, our simulation results demonstrate that the proposed FF 
approach provides substantial performance improvements as compared to the AF relay beamforming 
strategy. 
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II. Relay Network Model 



A. Filter- And- Forward Relaying Protocol 

Let us consider a half-duplex relay network with one single-antenna transmitting source, one single- 
antenna receiver (destination) node and R single-antenna relay nodes. Similar to [14] and [16]-[18], it is 
assumed that there is no direct link between the transmitter and destination nodes and that the network 
is perfectly synchronized. Each transmission is assumed to consist of two stages. In the first stage, the 
transmitting source broadcasts its data to the relays. The signals received at the relay nodes are then 
passed through the relay FIR filters to compensate for the effects of the transmitter-to-relay and relay- 
to-destination frequency selective channels. This type of relay processing corresponds to our proposed 
FF relaying protocol; see Fig. Q] As FIR filters have been commonly used for channel equalization in 
point-to-point communication systems, the FF strategy appears to be a very natural extension of the AF 
protocol to frequency selective relay channels. However, an important difference between these two cases 
is that in relay networks, it is meaningful to use a separate FIR filter at each relay node, while in the 
traditional point-to-point case, such a filter is commonly employed at the receiver. 

In the second transmission stage, the outputs of each relay filter are sent to the destination that is 
assumed to have the full instantaneous CSI. Using this knowledge, the receiver determines the filter 
weight coefficients of each relay according to a certain beamforming criterion. It is also assumed that 
there is a low-rate feedback link from the destination to each relay node that is used to inform the relays 
about their optimal weight coefficients. Alternatively, if the full instantaneous CSI is available at the relays 
rather than the destination, each relay node can determine its own weight coefficients independently^ In 
the latter case, no extra receiver-to-relay feedback is needed. Note that quite similar assumptions have 
been used in [14]-[18] in the frequency flat fading case. 

B. Signal Model 

Let us model the transmitter-to-relay and relay-to-destination channels as linear FIR filters 



'Note that network beamforming is commonly referred to as "distributed" because it is assumed that no relay can share its 
received signals with any other relay. 




(1) 
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where 

f; = [fi,i, ■ ■ ■ i fi,n] T ( 2 ) 

Si = [9i,i, ■ ■ ■ , 9i,r] T (3) 

are the R x 1 channel impulse response vectors corresponding to the Ith effective tap of the transmitter- 
to-relay and relay-to-destination channels, respectively. Here, f(w) and g(w) are the R x 1 vectors of 
channel frequency responses, and Lf and L g are the corresponding channel lengths, respectively. The 
Rx 1 vector r(n) = [ri(n), • • • ,r R (n)] T of the signals received by the relay nodes in the nth channel 
use can be modeled as 

Lf-l 

r(n) = ^ fis( n ~ + V{n) (4) 
1=0 

where s(n) is the signal transmitted by the source, rj(n) = [r)i (n), • • • ,rjji(n)] T is the R x 1 vector of 
relay noise, and (-) T denotes the transpose. Introducing the notations 

F — [foj • • • ,h f -i] 
s(n) = [s(n), s(n - 1), • • • , s(n — Lf + l)] T 

we can write (0]) as 

r(n) = Fs(n) + r}(n). (5) 
The signal vector t(n) = \t\{n), ■ ■ ■ , tn(n)] T sent from the relays to the destination can be expressed as 

L w -l 

t( n ) = ^ Wf r(n - I) (6) 



1=0 



where 



W ; = diag{u; Zi i, • • • ,w l R } 



is the diagonal matrix of the relay filter impulse responses corresponding to the Ith effective filter tap of 
each relay, L w is the length of the relay FIR filters, (•) denotes the Hermitian transpose, and for any 
vector x, the operator diag{«} forms the diagonal matrix containing the entries of x on its main diagonal. 
Correspondingly, for any square matrix X, the operator diag{X} forms a vector whose elements are 
the diagonal entries of X. Note that if L w = 1, then the FF transmission in (J6J) reduces to the AF one. 
Inserting (O into ©, we have 

L„-l 

t(n) = £ Wf (Fs(n - I) + r,(n -I)). (7) 
1=0 
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Let us define 



s(n) = [s(n), s(n — 1), • ■ ■ ,s(n- Lf - L w + 2)] T . 

It can be seen that the vector s(n — I) is a subvector of s(n). Using this observation, d7]) can be rewritten 

as 



L„-l 



t(n) = Wf(F / s(n)+r 7 (n-/)) 



z=o 



where 



l columns 



(L w — 1— I) columns 



= [Oflxl, ' • • j O-Rxl, F, O^xi, • • • , O^xi ], I = 0, • • • , L w — 1. 

Let us also define 



W 4 [W ,--- ,W L _!] T 



L w —1 columns 





F 




fo 


fx ... 


fi^-l OiJxl • 


• Orxi 




Fi 




OrxI 


fo fi 


••• f i/-l • 


■ Orxi 








OrxI 


Orxi ••• 


fo fi • 





IT 



(8) 



f)(n) ^ fo J (n), t; 7 (n — 1), ■ ■ ■ ,rj (n — L w + 1)] 
where OjvxM is the N x M matrix of zeros. Using these notations, we obtain that 

L„-l 
i=0 

and, therefore, ([8]) can be expressed as 

t(n) = W^^s(n) + W H fi(n). (9) 
The received signal at the destination can be written as 

Lg-l 

y{n) = gft(n-l) + v(n) (10) 

1=0 

where v(n) is the receiver noise waveform and is the channel impulse response vector defined in ([3]). 
Using (|9]), we can rewrite (fTOl) as 

Lg-l Lg-l 

y(n) = gfW^:Fs(n - I) + gfW H f)(n - I) + v{n). (11) 



1=0 



1=0 
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Taking into account that the matrices are all diagonal and using the properties of the Kronecker 
matrix product, we obtain that 

= [wfG,,..X-i G 'l 

= w H {I Lw G t ) (12) 

where 

A r T T lT 

wi = diag{Wj 
Gi = diag{a} 



Ijv is the N x N identity matrix, and (8) denotes the Kronecker product. Using (1121 . we can further 
rewrite (fTTT) as 

L g -1 L g -1 
y(n) = ^ w H (I Lw ® Gi)T~s{n -0+1] w H (I Lro &> G { )i/(n - + v(n). (13) 
z=o 1=0 

Defining 

s(n) = [s(n), s(n - 1), • • • , s(n - Lf - L w - L g + 3)] T 

fj(n) 4 [^(n), r/ T (?i - 1), • • • , r/ T (n - L w - L g + 2)] T 

we notice that s(n — Z) and r)(n — Z) are subvectors of s(n) and r/(n), respectively. Therefore, we can 
express ([T3l as 

L g -i L g -i 

y(n) = ^fc ® G)^zs(n) + ]T w h (I Lto 8) G,)Ii»K«) + w(n) (14) 
z=o Z=0 

where 

Z columns (L g — l—l) columns 



= [0RL ro xl, ■ • • , OflL ra xl, -T 7 , Orl„xi, ■ ■ ■ i Orl w xi], Z = 0, • • • , Lg — 1 
Z blocks (i g -l-Z) blocks 

I« = [Oi?L„,xK,- • ■, 0_rl iu xK, I_RL„, OijL raX _R,- ■ ■, OrL w xr], Z = 0, • • • , L g — 1. 

To express (TT4l in a more compact form, we further define 

£?4 [l£„®G 0j -" ,Il w ®Gl 9 -x] 

F = [^-,^ r i] T 
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and note that 

Lg-l 

^2(I Lw ®G l )F l = GF 

1=0 

L a -l 

l/.„ X G;)I/ : 

Z=0 

Using the latter two equations, (fl4l) can be expressed as 



Eft 



y (n) = w^^Fs(n) + w^C/Ir^n) + u(n) 



(15) 



w H £?[f,F] 



+ w H ^lT7(n) + u(n) 



Let f and F denote the first column and the residue of F, respectively, so that F = [f, F]. Then, (031 ) 
yields 

s(n) 
s(n) 

= w H g{s(n)+w H g~Fs(n)+\v H giri(n) + u(n) 
v v ' v s- 

signal ISI 

where 



(16) 



s(n) = [s(n - 1), • • • ,s(n - Lf - L w - L g + 3)] T . 

In (fl6l) . we can identify the three components 

1/5 (n) = w H Sfs(n) (17) 

yi (n) 4 w H 0Fs(n) (18) 

y n (n) = w^I^n) + t>(n) (19) 

as the destination signal, ISI, and noise components, respectively. Note that for the sake of computational 
simplicity of our techniques developed in the next section, block processing is not considered here, that 
is, the signal copies delayed by multipath are not coherently combined. 
The signal component in (fT71) can be expressed as 

y s (n) = w^G fos(n) 

= w^(g o 0f o )s(n) 

= w^hos(n) (20) 

where 

h o = g O 0fo (21) 
and denotes the Schur-Hadamard (elementwise) matrix product. 
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III. FlLTER-AND-FORWARD RELAY BEAMFORMING 

In this section, we develop three distributed FF beamforming approaches that utilize several alternative 
criteria. Our first FF beamforming technique is based on minimizing the total relay transmitted power 
subject to the destination QoS constraint, while our second and third approaches are based on maximizing 
the destination QoS subject to the total and individual relay transmitted power constraints, respectively. 
A useful modification of our third approach is also discussed, that enables to combine the later two types 
of constraints. 

A. Minimization of the Total Relay Power Under the QoS Constraint 

We first consider the distributed FF beamforming problem that obtains the relay filter weights by 
minimizing the total relay transmitted power P subject to the destination QoS constraint. As mentioned 
above, the destination QoS is given by the receiver SINR valudj and, therefore, the latter problem can 
be written as 

minP s.t. SINR > 7 (22) 

w 

where 7 is the minimal required SINR at the destination. 
Let us use the following two common assumptions 

E{s(n)s H (n)} = P s I Lf+La _ 1 , E{fi(n)f, H (n)} = tfl RLvj (23) 

on statistical independence of the signal and noise waveforms, respectively. Here, P s is the source 
transmitted power and is the relay noise variance. Using (O and (1231) . the transmitted power of 
the mth relay can be written as 

p m = E{|t m (n)| 2 } 

= E{eiW H Fs(n)s H (n)F H We m } + E{e^W ff jj(n)^(n)We m } 

= P s eJ n W H FF H We m + a 2 ri e T m W H We m (24) 

where e m is the mth column of the identity matrix. 

Using E m = diag{e m } and the properties of the Kronecker product, (1241 ) can be rewritten as 

p m = P s w H (I Lw E m ) TT H (I Lw 8) E m ) H w 

+a*w H E m ) (I £ „ (8) B m f w. (25) 

2 This is true because the processing at the destination is rather simple; in particular, no block processing is used. 
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The total relay transmitted power can be then expressed as 



where 



R / R \ 

P = Yl Pm = W " [Yl Dm W = W ^ DW (26) 
m=l \m=l / 



D m 4 p s (I Lm ® E m ) JFJF^ (I Ltu 8) F m f + cr 2 (I Lro 8) E m ) (I Lw ® E m ) H 
R i? 

D 4 £ D m = P s J] (1^ E m ) TT H (l Lw ® E m ) H + a 2 I flL „ . 



m=l m=l 

The SINR at the destination can be written as 



Using d20l) , we obtain that 



E{| yi (n)|2}+E{|y n (n)|2} 



E{|y s (n)| 2 } = E{|w "h s(n)| 2 } 
= P s w^h h^wo 
= P s w H A H h h^Aw 

= w H Q fl w (28) 



where 



A = \Ir>Qrx(l w -i)r] 
Q s 4 p s A H h h$A. 



Using (U8I ). we have 



E{\ yi (n)\ 2 } = E{w H gFs(n)s H (n)F H g H w} 



= P s w H gFF H G H w 

= w^Q.w (29) 

where 

Q, 4 P S QFF H Q H . 
Making use of ( fT9l , we also obtain that 

E{|y n ,(n)| 2 } = E{w^gif7(n)^(n)F^w} + C 7 2 

= w H Q n w + a 2 (30) 
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Q n ±a 2 v gu H g H . 



where 



Using (1261 ) and (I28l)- (l30b . the problem in (1221) can be rewritten in the following form: 

H n . W H Q S W 

mm w JJw s.t. — 77— 77— > 7. (31) 

w W H Q jW + Q„W + (T£ 

Introducing 

w^DVV Q^D^/ 2 (Q s - 7 Q i - 7 Q n )D~ 1 / 2 (32) 



we can reformulate the problem in (1311) as 

min ||w|| 2 s.t. w H Qw > 7a 2 . (33) 

w 

The constraint function in 03l ) can be used for checking the feasibility of the problem for any given 
value of 7. In particular, for all the values of 7 that lead to negative semidefinite Q, the problem in 
(133T ) is infeasible. It can be also easily proved that the constraint in (133T ) can be replaced by the equality 
constraint w^Qw = 7a 2 . Hence, the problem d33l ) is equivalent to 

min ||w|| 2 s.t. w ff Qw = 7a 2 . (34) 

w 

The solution of (l34b can be found by means of the Lagrange multiplier method. Let us minimize the 
Lagrangian 

H(w,X) = w H w + A(7cj 2 - w^Qw) (35) 



where A is a Lagrange multiplier. Taking gradient of (1351 ) and equating it to zero, we obtain that the 
solution is equal to that of the following eigenvalue problem: 

Qw = — w. (36) 
A 

Multiplying both sides of (l36l ) with Aw^ yields 

||w|| 2 = w H w = Aw^Qw = A7CT 2 . (37) 

It can be seen from (|3"7T ) that minimizing ||w|| 2 leads to the smallest positive A, which is equivalent to 
the largest 1/A in (l36l ). Using the latter fact, we conclude that the optimal solution to (I33T ) can be written 

as 

w opt = /3P{Q} (38) 
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where P{-} denotes the normalized principal eigenvector of a matrix and 

_/ 7fJ 2 y/a 
P ~ \nQ} H QHQ}) ' (39) 

Therefore, the optimal beamformer weight vector and the minimum total relay transmitted power can 
be expressed as 

w opt = 0I)-V 2 P{Q} (40) 

Pmin = 7^/Anax{Q} (41) 

respectively, where C max {-} denotes the largest (principal) eigenvalue of a matrix. 

Hence, the FF distributed beamforming problem (1221 enjoys a simple closed-form solution based on 
the principal eigenvector of the matrix Q. 

B. QoS Maximization Under the Total Relay Power Constraint 

Now, let us consider another useful distributed beamforming problem. Let us maximize the receiver 
SINR under the constraint that the total relay transmitted power does not exceed some maximal value 
P max . This problem can be written as 

max SINR s.t. P < P max . (42) 



max fln.,„x^p — s± w Dw - Pm - (43) 



Using d26l ) and d28l)-([30l), the latter problem can be expressed as 

w g Q s w ^ + mI1 

w W H QjW + w H Q n w + o% 

Introducing 

Q s ^ D- 1 /2 QsD -i/2 ) Q j+n A d-^CQ, + Q n )D~ 1/2 



we obtain that the problem (1431 ) can be rewritten as 

maX - rr~ Z ~ s -t- ll w ll < ^max (44) 

w w H Q i+n w + ai 

where, as before, w = D 1 / 2 w. It can be easily proved that the objective function in (l44l) achieves its 
maximum when the constraint is satisfied with equality (i.e., ||w|| 2 = P max )- Therefore, the problem ((44)) 
can be rewritten as 

w^Q s w ll2 
max = s.t. w = P max - (45) 

w W^(Q i+n + (a2/P max )I)w 

In contrast to the problem of Section ITlI-Al the problem (I43T ) is always feasible because for any positive 
Pmax, its feasible set is nonempty. Using the results of [21] (where a mathematically similar problem 
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has been discussed in a different context), we conclude that the objective function in (1451) is maximized 
when w is chosen as the normalized principal eigenvector of the matrix (Qi+ n + o"^/P max I) _1 Q s . Note 
here that any arbitrary scaling of w does not change the value of the objective function in ( 1451 ). However, 
the so-obtained vector w have to be properly scaled to satisfy the power constraint ||w|| 2 = P max - Then, 
the solution to (|45T ) can be written as 

w opt = y/K^P {(Q i+ „. + (vl/PmaxW^Qs) (46) 

and, therefore, the optimal beamforming weight vector and the maximum SINR at the destination can 
be written as 



w, 



'opt = V^D 4 ^ {(Qi+n + (vl/PmaxW^Qs} (47) 

SINR max = £ max {(Q l+n + (al/P^)!)- 1 ^} (48) 



respectively. 



C. QoS Maximization Under the Individual Relay Power Constraints 

Now, let us consider another relevant distributed beamforming problem which differs from (|42l is that 
the individual relay power constraints are used instead of the total relay power constraints. This problem 
can be written as 

max SINR s.t. p m < p m . max , m = 1, • • • , R (49) 

w 

where p mimax denotes the maximal transmitted power of the mth relay. Using d25l ) and (f28T>-([30T>. and 
introducing a new auxiliary variable r > [22], the problem (l49l can be transformed to 

max r 

W,T 

w ff Q,w 9 

w H D m w < p m , max , m = 1, • • • ,R. 

The first constraint in (l50l ) can be rewritten as 

^/F s \w H h\ > r^w H QiW + w H Q n w + erg (51) 

where h = A^ho- We observe that any arbitrary phase rotation of w will not change the value of the 
objective function in (l50l) . Using a proper phase rotation, we have that the constraint (I5TT) is equivalent 
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to 



^/F s w H h > r^/w^Q i w + w^Q ri w + a2 
Imjw^h} = 



(52) 
(53) 



where Im{-} denotes the imaginary part of a complex value. Note, however, that (|53l can be omitted as 
it is automatically taken into account in (l52l by the fact that the right-hand side of (l52l is non-negative. 
Then, the problem (l50l ) can be rewritten as 



max r 

w,r 



s.t. V^w^h > tJ w h QjW + w H Q n w + a% 
w H ~D m w < p m>max , m = !,-■■ ,R. 



(54) 



Let 



B 



T 



ORL„xl Qi + Qn 



V^V m , m = 1, 



,R 



be the Cholesky factorizations of the matrices B and D m , respectively. Introducing new notations 



w^[l,w T ] T , V m ^[0^„ xl ,V m ], h^[0,h 



TiT 



(55) 
(56) 

(57) 



we can further rewrite the problem d54l ) as 



max r 



s.t. V^w^h > r||Uw| 
||V m w|| < y/p 

m,maxi 



w H ei = 1. 



to = 1, 



(58) 



In contrast to the problem of Section MI- A I the problem (I58T ) is always feasible. This can be directly 
seen from its equivalent formulation d49l ) whose feasible set is always nonempty. Moreover, the problem 
(|58l ) is quasi-convex [22], because for any value of r, it reduces to the following second-order cone 
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programming (SOCP) feasibility problem: 

find w 

s.t. y^W s w H h > r||Uw|| 

||V m w|| < - V /Pm,max, TTl = 1, • • • , R (59) 
\v H e 1 = 1. 

Let r* be the optimal value of r in (I58T ). Then, for any r > r*, the problem (l59l is infeasible. On the 
contrary, if d59l ) is feasible, then we conclude that r < r*. Hence, the optimum r* and the optimal weight 
vector w* can be found using the bisection search technique discussed in [17]. Assuming that r* lies 
in the interval [77, t u ], the bisection search procedure to solve (l58l) can be summarized as the following 
sequence of steps: 

1) r := {n + t„)/2. 

2) Solve the convex feasibility problem d59"l ). If d59l is feasible, then 77 := r, otherwise r u := r. 

3) If (t u — 77) < e then stop. Otherwise, go to StepQ] 
Here, e is the error tolerance value in r. 

Note that the feasibility problem d59l is a standard SOCP problem, which can be efficiently solved 
using interior point methods [23] with the worst-case complexity of 0((RL w ) 3 ^). The initial interval 
for the bisection search can be selected as [77, t u ] = [0, y/ SINR max (P max ) ] , where SINR max (P max ) can 
be obtained from (|48l ) by choosing P max = J2m=i Pm,max- This particular choice is motivated by the 
fact that the optimal SINR of (02]) always upper bounds the optimal SINR of (l49l 

Remark: It is worth noting that the total power constraint can be easily added to d59l just as one more 
second-order cone constraint 

|| Vw|| < \APmax 

where V^V is the Cholesky factorization of D. This allows us to directly extend the approach of d59l ) 
to a practically important case when both the individual and the total power constraints have to be taken 
into account [16]. 

D. Relationships Between the Proposed Methods and Earlier Techniques in the Flat Fading Case 

Let us now explore the relationship between the proposed three methods and the techniques of [14], 
[16] and [17] in the specific case when all the channels are frequency flat and each relay filter is just a 
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complex coefficient (Lf = L g = L w = 1). In the latter case, the transmitted power of the mth relay in 
(l25l) can be simplified to 

Pm = P s w£ E m f f if E^wo + o*wf E m E£w 

= a2 f,m\ w 0,m\ 2 + 0"^ \wO,m\ 2 (60) 

where <? 2 rn — P s \fo,m\ 2 - Then, the total relay transmitted power can be expressed as 

R 

P = Pm = W^DqWq (61) 

m=l 

where 

D = diag{a 2 f l , • • • , aj R } + Ir. 
The received signal at the destination (fT3T ) can be simplified to 

y(n) = w^(f go)s(n) + w^(g r/(n)) + v(n) . (62) 

signal noise 

As the channel is frequency flat, there is no ISI term in d62~l) . Hence, the SINR reduces to SNR, and it 
can be written as 



w^Q n0 w + a% 



SNR= "o^»"" (63) 



where 



Q s0 4 p s (f g0 )(f g0 )" = p s h h£ 

QA 2j- ri |2 i |2t 

n0 = a v diag{\go t i\ , ••• , |ffo,ij| )• 

Introducing the variables < a m < 1 (m = 1, ■ ■ ■ , R) and using them to express the relay powers as 

Pm — Q^mPm, max (64) 

we obtain from (l60l and (l64l that 

Pm,max 



F0,m = «m 



tf + Pslfaml 2 ' 

In [14], it is proposed to compensate the phases caused by the transmitter-to-relay and relay-to-destination 
channels by a proper choice of the phase of each u>o,m- This gives 



wo,m = a m J H e J (65) 
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where 6 m = arg/o,m + arggo.m- Inserting (1651 ) into (1631 ). we have 



^Em=l a m | fo,m90,m I \Jl^+p'"\fa. m \ 2 1 



SNR = i 5 -- fs-g (66) 

°w Z^m=l ^+P s |/ , m | 2 



From (1661) . it can be seen that our problem d49l ) in the considered particular case is equivalent to 

I/O 



^if-TZ^m^ ^+P,|/o,m| 2 
s.t. <«!,•• • ,aR < 1 . (67) 



It can be readily verified that (1671) and the problem in [14] are identical. Moreover, the problem of 
[16] extends that of [14] to the case when both the individual and the total power constraints are used. 
Therefore, in the flat fading case where the AF strategy is used instead of the FF one, our approach 
of Section Ull-C I reduces to that of [14] and, with the additional total power constraint added to d59l ), it 
reduces to that of [16]. 

To understand the relationship of our three approaches of Section HIT] and the techniques of [17], we 
note that in the flat fading AF case the only difference between the problem formulations in [17] and 
our problem formulations is that an extra statistical expectation over all the random transmitter-to-relay 
and relay-to-destination channels has been used in [17]. Hence, in the flat fading case, the problem 
formulations of [17] transfer to our problems (l22l ). (l42l and (|49l ) when the instantaneous instead of the 
second-order CSI is used in the methods of [17] and when the FF strategy is replaced by the AF one in 
our techniques. 

IV. Simulation results 

In our simulations, we consider a relay network with R = 10 relays and quasi-static frequency selective 
transmitter-to-relay and relay-to-destination channels with the lengths Lf = L g = 5. The transmitter uses 
the binary phase shift keying (BPSK) modulation. The channel impulse response coefficients are modeled 
as zero-mean complex Gaussian random variables with an exponential power delay profile [24] 

= <y(t-zr.) (68) 

1=0 

where L x £ {Lf,L g }, T s is the symbol duration, <5(-) is the Dirac delta function, Pr is the average 
power of the multipath components, and at characterizes the delay spread. In our simulations, Pr = 1 
and at = 2T S are used. The relay and destination noises are assumed to have the same powers, and 
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the source transmitted power is 10 dB higher than the noise power. To obtain the bit error rate (BER) 
curves, it has been assumed that the symbol-by-symbol maximum likelihood (ML) decoder is used at 
the receiver. 

In the first example, we test the approach of (l40l which is based on minimizing the total transmit 
power subject to the QoS constraint. Fig. [2] displays the total relay transmitted power versus the minimum 
required SINR at the destination for different lengths of the relay filters. As for randomly generated signal, 
noise and channel values the feasibility of (l33l is a random event, this problem can be infeasible for some 
percentage of simulation runs. To deal with this fact in our first example, we call the problem ergodically 
infeasible when the number of simulation runs leading to infeasibility of (1331 is larger than the half of 
the total number of simulation runs; otherwise, this problem is classified as ergodically feasible. If the 
problem is ergodically infeasible, the corresponding points are dropped from the figures displaying the 
behaviour of the total transmitted power. In the case of ergodic feasibility, the corresponding points are 
computed by averaging over the "feasible" runs and displayed in these figures. For example, there are 
several dropped points in Fig. [2] at high SINR values that correspond to the case of ergodic infeasibility 
of d33l ). To further illustrate the effects of the required SINR and L w on the feasibility of the problem 
(l33T ). the probability that this problem is feasible is displayed in Fig. [3] The latter probability is referred 
to as feasibility probability. 

It can be seen from Fig. |2]that using the FF strategy at the relays, one can substantially reduce the total 
relay transmitted power as compared to the AF approach. Also, from Figs. [2] and [3] is is clear that the 
FF strategy significantly improves the SINR feasibility range of the considered distributed beamforming 
problem. These improvements are monotonic in L w : for example, according to Fig. |2j for SINR = 12 dB 
this problem is ergodically infeasible for the relay filter lengths L w = 1,2, but it becomes ergodically 
feasible for any L w > 3. These observations are further supported by Fig. [3] that demonstrates that the 
feasibility probability can be substantially improved by increasing the relay filter length L w . 

Figs. |4] and [5] depict the total relay transmitted power and feasibility probability versus the relay filter 
length L w for different values of the required SINR at the destination. Similarly to the previous two 
figures, Figs. |4] and [5] clearly demonstrate that the performance (in terms of the relay transmitted power) 
and feasibility of the QoS constraint can be substantially improved by using the FF approach in lieu of the 
AF strategy, and these improvements become more pronounced when increasing the relay filter length. 
Note that theoretically, to fully compensate the effect of frequency selective source-to-relay and relay- 
to-destination channels, it is required that L w > Lf + L g — 1. However, in the exponential power delay 
profile case (where these channels are mainly determined by several first taps), they are well compensated 
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even with L w < Lf + L g — 1. As follows from Fig. [4] depending on the SINR value, the filter lengths 
L w = 2 to L w = 5 appear to be sufficient. 

In the second example, we test the approach of d47l ) which maximizes the QoS subject to the total 
power constraint. Note that the problem (|45T ) is always feasible and, therefore, there is no infeasibility 
issue in this example. Fig. [6] shows the achieved SINR versus the total relay transmitted power P max for 
different lengths of the relay filters. Fig. [7] depicts the SINR versus the relay filter length L w for different 
values of the total relay transmitted power P mAx . It can be seen from these figures that the QoS can be 
substantially improved by increasing the filter length L w . 

To illustrate the receiver error probability performance of the FF relaying approach based on the 
particular example of the distributed beamformer (l47l) . in Figs. [8] and [9] we display the receiver BERs 
versus the total transmitted power -P max and the relay filter length L w , respectively. It can be observed 
from these two figures that increasing the filter length, we can substantially decrease the receiver error 
probability. 

In our last example, the approach of Section IIII-CI is tested, which maximizes the QoS under the 
individual relay power constraints using a combination of d59l ) and bisection search. As in the previous 
example, the underlying problem (I58T ) is always feasible and, therefore, there is no infeasibility issue 
here. It is assumed that all the relays have the same maximal allowed transmitted power p max . Fig. [10] 
displays the SINR versus p max for different lengths of the relay filters. Fig. [TT] shows the SINR versus 
the relay filter length L w for different values of p max . 

The conclusions following from Figs. [10] and [TT] are quite similar to that following from Figs. [6] and 
[7J As the individual power constraints are tighter than the total one, the SINRs achieved for any value 
of the total transmitted power in Figs. [TOl and [TT] are a few dB's lower than that achieved in Figs. [6] and 
[TJfor the same value of Rp ma , x . 

Summarizing, all our examples clearly verify that the proposed FF strategy substantially outperforms 
the AF approach in the frequency selective fading case. 

V. Conclusion 

The problem of distributed network beamforming has been addressed in the case when the transmitter- 
to-relay and relay-to-destination channels are frequency selective. To compensate for the effects of these 
channels, a novel filter-and-forward relay beamforming strategy has been proposed as an extension of the 
traditional amplify-and-forward protocol. According to the former strategy, FIR filters have to be used 
at the relay nodes to remove ISI. 
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Three relevant half-duplex filter-and-forward beamforming problems have been formulated and solved. 
Our first technique minimizes the total relay transmitted power subject to the destination QoS constraint, 
whereas the second and third methods maximize the destination QoS subject to the total and individual 
relay transmitted power constraints, respectively. For the first and second approaches, closed-form beam- 
formers have been obtained, whereas the third beamformer is computed using convex optimization, via 
a combination of bisection search and second-order cone programming. It has been also shown that the 
latter convex optimization-based relay beamformer can be easily extended to the practically important 
case when the individual and total power constraints should be jointly taken into account. 

Our simulation results demonstrate that in the frequency selective fading case, the proposed filter- 
and-forward beamforming strategy provides substantial performance improvements as compared to the 
commonly used amplify-and-forward relay beamforming approach. 
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Fig. 1. Filter-and-forward relay network. 
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Minimal required SINR (dB) 
Fig. 2. Total relay transmitted power versus required SINR; first example. 
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Minimal required SINR (dB) 



Fig. 3. Feasibility probability of the problem d33b versus required SINR; first example. 




Fig. 4. Total relay transmitted power versus relay filter length L w ; first example. 
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Fig. 5. Feasibility probability of the problem 03t versus relay filter length L m \ first example. 
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Fig. 6. SINR versus the maximal total relay transmitted power P max ; second example. 
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Fig. 7. SINR versus relay filter length L w ; second example. 
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Fig. 8. SER versus the maximal total relay transmitted power P max ; second example. 
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Fig. 9. SER versus relay filter length L w ; second example. 




Fig. 10. SINR versus the maximal individual relay transmitted power p max ; third example. 
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Fig. 11. SINR versus relay filter length L w \ third example. 
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